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In this paper we present a 5-parametric family of static asymptotically fiat solutions 
for the superposed gravitational and electromagnetic fields of two Reissner-Nordstrom 
sources with arbitrary parameters — masses, charges and separating distance. A pro- 
cedure for solving of the linear singular integral equation form of the electrovacuum 
Einstein - Maxwell equations for stationary axisymmetric fields is described in detail. 
$-H , The 4-parametric family of equilibrium configurations of two Reissner-Nordstrom sources 

bX)i (one of which should be a black hole and another one - a naked singularity) presented in 

our recent paper^ arises after a restriction of the parameters of the 5-parametric solution 
presented here by the equilibrium condition which provides the absence in the solution 
^ of conical points on the symmetry axis between the sources. 

Introduction 

' In our recent short paper^ we presented (in a surprisingly simple form) an exact 4- 

parametric family of static asymptotically flat solutions of electrovacuum Einstein 
- Maxwell equations which describes the equilibrium configurations of two (non- 
. rotating) charged massess in General Relativity. In this paper we present a more 

general, 5-parametric solution of these equations which represent a nonlinear super- 
position of fields of two Reissner-Nordstrom sources with arbitrary mass and charge 
parameters and arbitrarily chosen separating distance and describe the procedure 
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5h ' for a construction of a superposition of these fields. In the subsequent sections, we 

^ use a special divergent form of the reduced Einstein - Maxwell equations to derive 

the Komar-like integrals for the total gravitational mass and charge of this field 
and calculate the physical masses and charges of the sources defined as the additive 
inputs of each source in the total gravitational mass and total charge of the system. 
The expression of our solution in terms of these physical parameters simplifies it 
considerably. Then we determine the constraint which should be imposed on the 
parameters of this superposition of fields to provide the absence in the remaining 
4-parametric solution of any non-physical singularities such as the conical points on 
the axis outside the sources. This constraint plays the role of the condition for equi- 
librium of these sources in their common gravitational and electromagnetic fields. 

The 5-parametric solution, presented in this paper, was constructed using the 
monodromy transform approach'^'^ which give rise to a reformulation of the Einstein 
- Maxwell equations for stationary axisymmetric fields in terms of equivalent system 
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of linear singular integral equations. As it is explained below, the functional param- 
eters (monodromy data) in the kernels of these integral equations for our solution 
are chosen as analytically matched, rational functions of the spectral parameter. 
This choice of the monodromy data means that the solution possess rational axis 
data, i.e. that the values of all components of metric and electromagnetic potential, 
expressed in terms of cylindrical Weyl coordinates {p,z), on the axis of simmetry 
p = are rational functions of z. This means also that our solution virtually is the 
soliton solution on the Minkowski background in the sense, that it consists of a pure 
soliton part and its analytical continuation in the space of parameters^ 



On the integral equation methods 

Different approaches to solution of integrable space-time symmetry reductions of the 
Einstein's field equations, which began to develop about thirty years ago, gave rise 
to different formulations (generally non-equivalent to each other) of the integral 
equation methods, which allow to calculate the solutions of these field equations 
solving some systems of linear singular integral equations. 

The first reformulation of symmetry reduced Einstein equations in terms of 
a system of linear singular integral equations for vacuum gravitational fields (to- 
gether with a construction of vacuum solitons) was proposed in the framework of 
the inverse scattering approach'^. This construction was based on a formulation of 
equivalent Riemann - Hilbert problem for the 2 x 2-matrix functions on the spec- 
tral plane which gave rise to some system of linear singular integral equations. It 
is important, that the solutions of this equations possess the character of solution 
generating transformations, because the matrix kernel of these integral equations 
includes, as functional parameters, the components of an arbitrarily chosen vacuum 
metric which serves as the background for solitons, while the generating solution it- 
self can be considered as describing some nonlinear perturbation of this background. 

In the framework of another, group-theoretic approach, for construction of the 
solution generating transformations corresponding to the elements of the infinite di- 
mensional algebra of internal symmetries of electrovacuum Einstein-Maxwell equa- 
tions for stationary axisymmetric fields (found by Kinnersley and Chitre^), Hauser 
and Ernst^ reduced these equations to a homogeneous Hilbert problem for 3x3- 
matrix functions of an auxiliary complex parameter and then, to the corresponding 
3 X 3- matrix linear singular integral equations. This construction assumed an ad- 
ditional constraint imposed on the class of solutions, which means that only those 
solutions are considered which possess (locally) a regular behaviour of fields near 
the axis of symmetry. The integral equations^ have rather complicate matrix kernel 
which construction includes a calculation of matrix exponents of the elements of 
Kinnersley and Chitre algebra, represented by some algebraically defined holomor- 
phic 3 X 3-matrix functions of the mentioned above auxiliary complex parameter. 



*This analytical continuation is similar to that which connects different parts of the Kerr-Newman 
family of solutions corresponding to the field of a naked singularity and the black hole solution. 
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Later, fixing the choice of the seed solution by the simplest one, SibgatuUin^ reduced 
the matrix integral equations® of Hauser and Ernst to a simpler scalar linear singu- 
lar integral equation with some "normalization" condition imposed additionally on 
its solutions. The kernel of this scalar integral equation was expressed explicitly in 
terms of the values of the Ernst potentials on the axis of symmetry. This integral 
equation was actively used during the last one and a half decades by SibgatuUin 
and others for mostly formal calculation of asymptotically flat solutions for various 
particular choices of (asymptotically flat) rational axis data for the Ernst potentials. 

The monodromy transform approach^'^ does not follow the ideology of the matrix 
Riemann-Hilbert problems, however, it is also based on some ideas of the modern 
theory of integrable systems, analytical theory of differential equations and the the- 
ory of linear singular integral equations. For physically different classes of vacuum 
and electrovacuum fields with two commuting isometrics (stationary axisymmetric 
fields, plane and cylindrical waves, inhomogeneous cosmological solutions and some 
others) this approach suggests rather simple general construction of the coordinate- 
independent functional parameters (called as monodromy data) which characterize 
uniquely every local solution. The problem of constructing solutions for given mon- 
odromy data (the inverse problem of the monodromy transform) gave rise to a 
system of linear singular integral equations which differs essentially from the in- 
tegral equations mentioned above. A specific structures of the (scalar) kernels of 
these integral equations and of the integration path on the spectral plane allow 
to describe all degrees of freedom of the gravitational and electromagnetic fields 
and make these integral equations equivalent to the symmetry reduced Einstein - 
Maxwell equations. For stationary axisymmetric fields satisfying the regularity axis 
condition, the integral equations^''^ simplify considerably. A general scheme for con- 
structing of solutions with any rational analytically m,a,tched m,onodromy data (or, 
equivalently, of solutions for any rational, not only asymptotically flat, axis data) 
was described earlier in detail (see^'^ and the references therein). 

Solitons on the Minkowski background and rational axis data. 

A discovery of existence of pure gravitational solitons and of the ways for their 
generating on arbitrarily chosen vacuum backgrounds - the dressing methods'* (see 
also^), together with initiated by these results developments of the similar methods 
for Einstein - Maxwell fields*^, gave us a powerful tool for construction of a large 
variety of solutions of Einstein's field equations and for nonlinear superposition of 
fields of certain kinds of sources with various external fields. 

On the other hand, later developments of the integral equation methods briefiy 
described above suggested a simple idea to calculate the solutions of Einstein and 
Einstein - Maxwell equations using these integral equations with the simplest choice 
of functional parameters (the contour data, or the axis data for the Ernst potentials, 
or analj^ically matched monodromy data) in their kernels as rational functions of 
their arguments. This also leads to a construction of large families of exact solutions 
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with arbitrary (finitely large) number of free parameters, but without any freedom 
in the choice of the background solution. This gives rise to the obvious questions 
concerning a comparison of these two constructions. 

First of all, it is clear, that the asymptotically flat rational-data solutions can be 
only a very special case of soliton solutions corresponding to a particular choice of 
the background for solitons. This can be the Minkowski space-time or any other soli- 
ton solution on this background. Indeed, a direct calculation shows, that vacuum 
solitons''^ and electrovacuum solitons^°, both generated on the Minkowski back- 
ground, represent the asymptotically flat solutions with rational axis data. These 
solitons have as many of free parameters per one simple pole of the dressing matrix, 
as it is necessary to describe an arbitrary chosen rational asymptotically flat axis 
data with a twice lower number (for vacuum case) or with the same number of poles 
in this dataQ Of course, the coincidence of the number of parameters itself does not 
mean that these classes of solutions coincide. Moreover, some doubts in this coinci- 
dence can arise from that fact that the construction^" of electrovacuum solitons (in 
contrast to vacuum soliton generating technique"*) allows to generate the solitons 
with arbitrarily located complex poles of the dressing matrix and with the complex 
conjugated poles of the inverse matrix, while the solitons with real poles do not 
arise in this technique. Of course, one can try to obtain the electrovacuum solitons 
with real poles as limiting cases of solitons with complex poles, however, this leads 
to solutions with coinciding real poles of dressing matrix and its inverse, which have 
less number of free parameters per each pole than the solitons with complex poles. 

To clarify this situation, it is worth to consider a one-soliton solution on the 
Minkowski background, which coincides with the over-extreme part of the well 
known Kerr-Newman solution corresponding to a naked singularity. The under- 
extreme part of this solution, which has a horizon and corresponds to a black hole, 
can be described in terms of the soliton generating technique by a dressing matrix 
with one real pole. However, as it is well known, the under-extreme part of the 
Kerr-Newman family can be obtained by a simple analytical continuation from its 
over-extreme part in the space of its parameters. It is very likely, that the sim- 
ilar analytical continuation is possible not only for a one-soliton solution on the 
Minkowski background, but for any number of solitons as well. Indeed, solving the 
integral equations for the (static) two-pole monodromy data, we have found that 
the part of our static 5-parametric solution, which corresponds to a pair of naked 
Reissner-Nordstrom singularities, coincides with the static subfamily of two-soliton 
solutions on the Minkowski background, while the remaining part of this family 
is connected with this two-soliton solution by the similar analytical continuations 
in the space of parameters as one uses in the one-soliton case to connect different 
parts of the Kerr-Newman family of solutions. The motivation given above shows 



^It is necessary to note, however, that the poles of the dressing matrices for solitons coincide with 
the poles of the solutions of the integral equations, but they do not coincide in general, with the 
poles of the axis data for the Ernst potentials. 
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that there is the most reason to expect that the class of solutions of the integral 
equations for asymptotically flat rational axis data consists of the already known 
solitons on the Minkowski background together with their analytical continuations 
in the space of parameters (which should be supplied also, for completeness, by cer- 
tain limiting solutions corresponding to multiple poles in the dressing matrices) . On 
the one hand, in view of the above considerations, this means that there is no any 
reason to consider every formally calculated solution of the integral equations for 
rational asymptotically flat axis data as some new one describing the "extended" 
or "generalized" solitons, because it would be curiously enough to call similarly 
the under-extreme part of the Kerr-Newman family of solutions which is a simple 
analytically continuation in the parameter space of its over-extreme part. 

On the other hand, the above motivation does not mean at all, that there is 
no sense to use the integral equation methods for calculation of these soliton solu- 
tions starting directly from some asymptotically flat rational axis data. It is worth 
to mention here, that the application of the soliton generating technique leads to 
some specific parameterization of the constructing solutions which differs signifi- 
cantly from that which arise for the solutions constructed from the corresponding 
solutions of the integral equations. Each of these parameterizations can occur to be 
more or less useful in different considerations of this class of soliton (or rational- 
data) solutions. For example, it seems, that just the existence of the analytical 
continuation in the parameter space discussed above can be proved for N-soliton 
solution more easily if we use its expression which arises from the solution of the in- 
tegral equations, but considering the general twelve-parametric two-soliton solution, 
we find that it has the most compact form just in the soliton parameterization.^ 

This concludes our very brief reminding of some fragments of the history of 
the methods and of some interrelations between the solitons and solutions with 
rational axis data derived from the integral equation methods. At the end of this 
Introduction, we sketch out the monodromy transform approach and the integral 
equation method which we use for construction of our 5-paramctric sohition for a 
superposition of fields of two Reissner-Nordstrom sources. In all necessary details, 
a theory of this method and a general scheme for construction of solutions with 
rational, analytically matched monodromy data, have been developed long ago and 
described in the papers cited above with some additional useful references therein. 

Monodromy transform approach 

This approach is based a) on the parametrization of the space of local solutions of 
the symmetry reduced Einstein - Maxwell equations by the monodromy data - a 
set of coordinate independent functions of a spectral parameter, which determine 
the branching properties on the spectral plane of the fundamental solution of as- 
sociated linear system, and b) on the reformulation of these equations in terms of 
a system of linear singular integral equations. For stationary axisymmetric fields, 
this approach allows to construct nonlinear superpositions of electrovacuum fields of 
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different sources cliaracterizcd by analytically matched, rational monodromy data^ 
(see also the Appendix in*). In the following sections, we demonstrate at first that 
the external field of a single Reissner-Nordstrom source is characterized just by this 
kind of the monodromy data functions which have on the spectral plane one simple 
pole and vanish at infinity. For superposition of fields of two such sources we choose 
the rational monodromy data functions as the sums of two poles with such coeffi- 
cients which guarantee that the solution is static. Given this monodromy data, we 
describe step by step the construction of the corresponding solution. 

The space of local solutions and its parameterization by 
monodromy data 

For electrovacuum Einstein - Maxwell fields depending only on two space-time coor- 
dinates, in the entire space of local solutions, which are analytical near some initial 
point and take at this point (the point of " normalization" ) some "standard" values, 
every local solution with the metric gik and electromagnetic potential Ai can be 
characterized uniquely by the monodromy data which consist of four coordinate- 
independent holomorphic functions of the spectral parameter w. 

{gik{x\x^),A,{x\x^)} < — > {u±(w),v±(w)} (1) 

This monodromy data are defined as a complete set of independent functions which 
characterize the branching properties of the corresponding fundamental solution of 
the associated linear system on the spectral plane at its four singular points. This 
data arc defined uniquely for any local solution and for arbitrary choice of these 
functions there always exists a unique local solution of electrovacuum Einstein - 
Maxwell equations with given monodromy data. 

This construction is useful as far as some effective methods can be developed for 
explicit constructing of local solutions for given monodromy data. Then, given an 
explicit local solution, the corresponding global solution can be determined using its 
analytical continuation to other space-time regions where this solution may reveal 
also some regular behaviour or approaches various types of singularities. 

The m,aster system, of linear singular integral eguations 

The key point of the use of the monodromy data is the existence of some system of 
linear singular integral equations whose kernels and right hand sides are expressed 
algebraically in terms of the monodromy data and whose solution determines (by 
means of some quadratures) all components of metric and electromagnetic potential. 

In general, the structure of this system of linear singular integral equations is 
rather complicate. The singular integrals are defined on the contour L on the spectral 
plane which consists of two disconnected parts L = _L+ + L_. The locations of their 
endpoints depend on the space-time coordinates and coordinates of the initial point 
which enter also the integrands as parameters. In particular, goes from w = 
to w = ^ and L_ goes from w = rio to w = r], where, for example, for stationary 
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axisymmetric fields in terms of Weyl coordinates £, = z + ip, rj = z — ip, but for 
plane waves £, = x + t,r] = x — t) and {£^0, Vo) correspond to the initial point. 

The monodromy data as well as the unknown functions in these integral equa- 
tions arc defined in two disconnected regions of the spectral plane the neighbour- 
hoods of L-)- and L_, where they are represented by pairs of functions (u+, u_), 
(v+, v_) and (cp+, cp_) respectively. It is clear that these equations can not be 
solved for arbitrarily chosen monodromy data. However, for some classes of fields, 
such as, for example, stationary axisymmetric fields with a regular axis of symmetry 
considered in this paper, these integral equations can be simplified considerably and 
admit infinite hierarchies of multiparametric families of explicit solutions. 

Monodromy data for stationary fields with a regular symmetry axis 

For stationary axisymmetric fields, it is typical that physical and geometrical for- 
mulations of various problems (as, for example, in the case of asymptotically flat 
fields) imply that at least some part of the axis of symmetry is free of the field 
sources and therefore, a behaviour of metric and matter fields in the neighbourhood 
of this part of the axis should be regular. For such fields, if we choose the initial 
point (the point of normalization) of a solution on such regular part of the axis, the 
initial points of the contours L+ and L_ coincide and instead of two disconnected 
contours we obtain one simple curve L = L+ + Changing preliminary the di- 
rection of integration on L_, we obtain one contour where the system of singular 
integral equations is defined. It starts at the point w = r], goes through the initial 
point IV = zo and ends at w = ^. On this contour the monodromy data and the 
unknown variable in the integral equations are represented by twice lower number 
of holomorphic functions because for these fields we have 

u+{iu) = U-{w) = u{iu), v+{w) = V-{w) = v{w), (p +{w) = -{w) = (p {w) 

(2) 

Usually, we call these conditions as the regularity axis condition, however we note 
that these conditions guarantee only a regular local behaviour of fields near the 
axis, but they do not exclude the presence of some non-curvature singularities, such 
as conical points on the axis or closed time-like curves near it. Thus, for stationary 
axisymmetric electrovacuum fields with a regular axis of symmetry the monodromy 
data are represented only by two arbitrary holomorphic functions u{w) and v(i(;). 
We recall also that v{w) is "responsible" for a presence of electromagnetic field, 
so that for vacuum v{w) = and the space of solutions of vacuum stationary 
axisymmetric fields near the regular part of the axis of symmetry is parameterized 
by the monodromy data consisting of one holomorphic function u{w). 

Another important property of the stationary axisymmetric electrovacuum fields 
with the regular axis of symmetry is that any solution can be characterized by the 
finite values of its metric components and potentials (e.g., of the complex Ernst 
potentials £{p, z), $(/), z)) on the regular part of the axis of symmetry. In this case. 
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the monodromy data ^ can be related to the values £{z), $(2) on the axis: 

£{p = 0, z) ^ £o-2i{z- zq)u{w = z), ^{p ^ 0, z) ^ ^o + 2i{z~ zo)v{w = z), (3) 

where z is the Weyl coordinate along the axis, zq is a coordinate of the initial point 
on this axis, £q and $0 are the "normalized" values of the Ernst potentials at the 
initial point, for which we usually put Sq = I and $0 = 0- 



The conditions for asymptotically flat and static fields 

The expressions ([3]) allow us to relate the structure of the monodromy data ([2]) with 
some physical and geometrical properties of fields. In particular, for any asymptot- 
ically flat field u{w) and ^{w) should be holomorphic at w — 00 and 

u(w) — !■ and v(u;) — s- for w 00. 

Moreover, in this case, the coefficients of expansions of these functions in the in- 
verse powers of w for w 00 can be simply related to the multipole moments of 
this asymptotically flat field. Therefore, the multipole structure of the field can be 
determined in advance by the appropriate choice of the monodromy data. 

For static fields, £ should be real, while $ should be real for pure electric fields 
and imaginary for pure magnetic fields and therefore, u(w) and v(w) should satisfy 

u^{w) — — u(w) and v^(w) = =Fv(i(;) 

where u^(w) = u{w) and v'^{w) = ^{w) and a bar means a complex conjugation. 



Exact solutions with rational monodromy data 

An infinite hierarchies of solutions of the Einstein - Maxwell equations can be cal- 
culated explicitly if we choose the analytically matched monodromy data ^ to be 
rational functions of the spectral parameter: 

U{w) V{w) 



Q{w)' ' ' Q{w) 

where the functions U{w), V{w) and Q{w) are some polynomials. A general algo- 
rithm for solution of the integral equations for these polynomials of arbitrary orders 
was described in.'^'^ This algorithm leads to explicit form of solutions in a unified, 
but rather complicate form, and, as we shall see below, a large careful work is nec- 
essary for finding of appropriate choice of physical parameters which can simplify 
significantly the constructed solutions. 



The linear singular integral equations 

For stationary axisymmetric fields outside their sources the metric and electromag- 
netic vector potential can be considered in cylindrical coordinates in the form 

ds^ = gttdt^ + 2gt^dtdip + g^^dtfP- - /(dp^ + dz'^) 

(4) 

A, = {At, 0,0, A^}, 
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where = {t,p,z,ip\ and the metric components gtt, Qtip, gtptp and / as well as 
the components A-t and of the electromagnetic potential are functions of the 
coordinates p and z only. It is well known, that for the metric and electromagnetic 
fields ([3]) the symmetry reduced electrovacuum Einstein - Maxwell field equations 
decouple into two parts. One of these parts is a closed system of "dynamical" 
equations - the nonlinear partial differential equations for the functions gu, Qtipi 
g>fi>fi, At and A^p. (We recall here that the Weyl cylindrical coordinates are defined 
so that gttgipip ~gtip = and therefore, only two of these three metric components 
are unknown functions) . Another part is a pair of constraint equations which allow 
to determine in quadratures the conformal factor / provided the solution for gtt-, 
gtip, gtptpi At and A^p is already known. Though an explicit calculation of these 
quadratures for the conformal factor / represent usually a large technical difficulty, 
the principal problem is a construction of the solution of the " dynamical" equations 
with wanted physical and geometrical properties. That is why we concentrate below 
mainly on the construction of solution of the "dynamical" equations. 

In accordance with a general scheme, ^'^ the components of metric and electro- 
magnetic potential (|4]), which satisfy the electrovacuum Einstein - Maxwell equa- 
tions, can be expressed in the form 

gtt = l- i{Rt'^ - Rt^) + *t^t £ = 1 - 2iRt'^ 

gtv = -iiz-za)+i{Rt* + R^f) + <^;^^ (5) 

g^^ = i{R^' - R^*) + $^1^ K^v) " ' UvV 

where {$t,$^} are the components of a complex electromagnetic potential and 
Re$t = At, Re$^ = A^; the functions Rl, Rf, Rl, i?^, i?^, i?* constitute a matrix 
which is determined by the integral over the contour L on the spectral plane 

_ fRt* Rf R*t\^]_ /l + *(C-^o)u^(C)^ 
i?^ R*J~i7:r \-t{(:-zo) 

(6) 

Here C, & L and the contour L (unlike the general case) is a simple curve which 
starts from w ~ rj = z — ip, goes through the initial point w — zq and ends at 
the point w — = z + ip, where p and z are the well known cylindrical Weyl 
coordinates in which gitg^pip — gtip = ^P^, and zq determines the location of the 
point of normalization on the axis of symmetry p = 0. In the above integral, [A]^ 
denotes a jump (i.e. a half of the difference between left and right limits) at the 
point C, £ L oi a, "standard" branching fmiction 



R^k%^;^:)=-/[Ak( _ / )®{<p['i(c),cpM(c),cpM(c)}dc 



\ = ^{C-z~ip){C~z + v)/(C - zo)2, A(p,z,C = (X3) = l. (7) 

The functions cp I-*"! (C), cp (C), tp (C) should satisfy to the decoupled linear singu- 
lar integral equations with the same scalar kernels and different right hand sides:^''^ 

>''i(C)\ / 1 \ 

(pM(C) dC= u(r) (8) 

,(pM(c)y Vv(r)y 



[A]c^(t,C) 



L 
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where (^,t Cz L; there is a Cauchy principal value integral in the left hand side, and 
the kernel function TI(t, C) in its integrand is 

H{t, C) = 1 + *(C - zo) [ut(C) - u(r)] + 4(C - zo)'v(r)vt (C) 

and everywhere below, we can put zq — without any loss of generality. The 
coordinates p and z enter the equations ([5]) as the parameters which determine the 
location of the endpoints of the contour L and as the arguments of the function A, 
however, for simplicity we have not shown explicitly in ([5]) the dependence of the 
unknown functions (p t"""' , cp I"! and cp ["^1 on these coordinates. We note also, that in 
general, (p , cp ["1 and cp ["^1 as well as the right hand sides of ([5]) constitute the 
row-vectors, however, for a convenience we write (|H1) in a transposed form. 



The Kerr - Newman field as a one-pole solution 

We begin our description of solution of the integral equation ([5]) with a simple case 
of a one-pole structure of the monodromy data functions: 

uM = r, V w) = -. (9) 

w — h w — h 

where Mq, vq and h are arbitrary complex constants. For these data the kernel 
function 7i(T, C)/(C — r) can be split into the singular and regular parts: 



H(r,C) 1 



(10) 



C-T {r-h){C-h) 

where the polynomials P(C) and i?(C) possess the expressions (zq — 0): 

P(C) = [1 - i{ua -u^) + 4woW] - C[(l + + (1 - iuo)h] + hh 

RiO = -il + iu^)C + h ^^^^ 

Assuming that h 0, what means that the pole is not located on the integration 
path L, we present the polynomial P{w) in a factorized form 

P{w) = Po{w - wi){w - wi), Po = 1 - i(tio -lio) +4woW (12) 

The coefhcients of P{'w) always are real and therefore, its roots wi and wi are real 
or complex conjugated to each other. We parameterize these roots as 

Wi — zi + ai, Wi — Zi — Gi (13) 

where zi is a real parameter while ai can be real or pure imaginary. A comparison 
of PT|) and allows to express zi and cti in terms of mq, Vq and h, however, it 
is more convenient to use zi and ai as new parameters and express some of the 
parameters Uqj ''^o ^^id h as functions of zi, ai and others. We give the explicit 
expressions later, but now we concentrate on solving of the integral equations ([5]). 

As it can be concluded from the structure of the singular integral equations, 
their solutions for the monodromy data ^ should have the form 

{cpW(u;), cpH(,„), (pM(uO} = ^^(^{^ + ^0, Y,, Zo} (14) 
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where Xq, Yq and Zq are independent of the spectral parameter w, but they can 
depend on the coordinates p and z. As we see from these expressions, the solutions 
of the integral equations ([H]) are rational functions of the spectral parameter and 
they have the poles coinciding with the roots of the polynomial P{w). 

To calculate the values of Xo, Yq and Zq explicitly, we substitute ([9]), pO]) and 
(fT4|) into the integral equations ([8]) and obtain the linear algebraic equations for 
Xq, Yq and Zo with rather complicate coefRcients. These coefficients are linear 
combinations of the usual or singular Cauchy principal value integrals of the form 

L L 

where t,(^ G L and fc > is some integer. Any integral of these types can be 
calculated explicitly using the elementary theory of residues. Indeed, the integrands 
in these integrals can be expressed as the jumps of analytical functions A(w) w'' /{w— 
t) and X{w) w'' / P{w) respectively. Therefore, these integrals can be expressed as 
the integrals over the closed curves C surrounding the contour L: 

ntj C-r 2mJ x-r mj P(C) ''^ 2m J P{x) ^ 

L C L C 

(15) 

where x ^ ^ ^'^d C surrounds L in negative direction (i.e. so that the interior is to 
the right) and close enough to L so that no poles of P{'w) are inside C. The function 
\{w) w'^/ {w — t) is analytical outside the contour L and it can have only the pole 
at w = oo, while the function X{w) w'^ / P{w) is analytical outside L, besides the 
poles which arise from zeros of P{w) and possible poles at w = oo. This allows us to 
transform the closed integration path C into a one approaching w — oo taking into 
account the inputs from the finite poles and therefore, in (fT5|) . each of the integrals 
over C is equal to the sum of residues of its integrand at finite poles plus the residue 
at w = CX3. Thus, for the integrals ((T5|) we obtain 

-/^rfc--E(A)fc^™r", 

1 f [A]cC^ ^ Xiwi)w'l Xiwi)w'[ _ f Xjw) 
mj P(C) ^ P'iwi) + P'iwi) [ P{w) 



(16) 



where (A)„i-fc means the coefficient in front of l/tu'^^™ in the inverse powers ex- 
pansion of the function X{w) at w — oo and similarly, (. . .)_i means the coefficient 
in front of l/w in the inverse powers expansion of a function at w = oo. The pa- 
rameters wi and wi are the roots of the polynomial P{w), while P'{wi) and P'{wi) 
mean a derivative of P(w) with respect to w at w = wi or w ^ wi respectively. 

The substitution of (fT4|) and (fTO|) into ^ and the subsequent calculation of 
the contour integrals leads to the equations whose left and right hand sides, after 
multiplication by {r—h), become polynomial functions of r. Equating the coefficients 
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of these polynomials, we obtain the linear algebraic equations for Xq, Yq and Zq 
whose solution allows to obtain the explicit expressions for cp : 

r [11. N [ul. N [vl/ n\ h + Ai uo vo\ , . 

jcpHH, cpl'H, cpHH| = -^|^ + z-^— , -| (17) 

where Aq and Ai possess the expressions 

^1 = + - + ^) (18) 

1 + ZMq r-p, ,-| h 

-2 — I - iuq) + h(l + iuo) \ + — . 

and the polynomials P{() and R{C) were defined in pT|) . Now we have to substitute 
these solutions of the integral equations into the integrands of quadratures ^ and 
calculate these quadratures using the same methods for calculation of these type of 
integrals as it was described above. The components of metric and of electromag- 
netic potential, as well as the Ernst potentials for the constructing solution can be 
calculated then algebraically from the expressions ([5]). 

The formal calculations for the monodromy data ([9]) described above lead us to 
the explicit but very complicate form of the solution. This solution can be simplified 
considerably using (i) a new and most appropriate set of parameters, (ii) some global 
gauge transformations - SL{2, i?)-rotations and rescalings of the Killing vectors 
and with the corresponding linear transformations of t and tp coordinates and 
(iii) a set of more convenient coordinates (such as the prolate or oblate ellipsoidal 
coordinates) instead of the Weyl coordinates p and z. In the remaining part of 
this section we describe these simplifications and show that this solution represents 
nothing more but the well known Kerr - Newman family. 

At first, comparing the explicit expression (fTTj) for P{w) with its factorized form 
(fT2|) and using (fT3|) . we obtain a set of relations between zi, a\ and the original set 
of independent parameters uqj ^^o a-nd h. Solving these relations we express uq) 'I'o 
and h in terms of a new set of parameters which consists of the real parameters m, 
a, 6, z\ and a complex parameter e: 



ih 

itQ = — i H —{zi + ia), h — zi + m — i{a + b) 



WiWi 



le 



hh 



(19) 







Zh WlWl 

and the parameter cti is is related to the parameters m, a, b, zi and e by the equation 

2 2 , i2 2 - 

ai — m +0 — a — ee 



Now, for the choice of new coordinates we note, that besides a presence of p and z 
in the solution explicitly, its dependence of these coordinates also comes from the 



February 2, 2008 3:44 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



13 



functions A(/o, z,k7) which should be calculated at the points w = wi and w = Wi, 
i.e. at the roots of P{w). Let us express these functions in the form 

Zl + (71 Zl — (Tl 

where xi and yi will be considered as new coordinates. If we compare these expres- 
sions with the definition ([7|), we obtain the relations between p, z and xi, yi. 



P= yxl-ol^l-yl, z^zi+xiyi 
Finally, we have to make some gauge transformation of our solution: 

1 H to = 1 - 2zuo 

$0 = 2iwo 



After this reparametrization, change of coordinates and gauge transformations we 
obtain the solution in a simple form. In particular, we obtain the Ernst potentials 

c ^ I 2 + 4eot;ot'j)) _ ^ 2{m - ib) 



htt - 


^ htt 
' Po 


£ 


htip - 


ht^p 











PqAq Xi + iayi + m — ib 

^ 2ivQ e 



/^Ao xi + iayi + m - ib 

which can be identified with those for the well known Kerr-Newman solution. The 
coordinates xi and yi are connected directly with the polar spheroidal coordinates: 

Xi = ri — m, yi — cos^i 

and the real parameters m, a, b, zi, real and imaginary part of e can be identi- 
fied respectively with the mass, angular momentum, NUT-parameter, electric and 
magnetic charges of the Kerr-Newman source (a black hole or a naked singularity). 



The Reissner - Nordstrom field as a static one-pole solution 

As it was already mentioned above, to obtain a static solution with the monodromy 
functions having only one pole, it is not necessary to make all of the calculations 
described above, but it would be enough to restrict from the beginning the mon- 
odromy data functions by the case, when these functions would take pure imaginary 
values on the real axis on the spectral plane. For this, it is necessary to choose the 
parameter h to be real and the parameters uq and vq pure imaginary. In this case, 
as one can see from (fT9| . we would have 

h — h, Uq = —Uq, Vq = —Vq <==^ = 0, 6=0, 6 = 6 

i.e. we obtain the Reissner - Nordstrom solution. 
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Superposing the fields of two Reissner - Nordstrom sources 

In this section we describe tlic key points of the calculation of solutions for mon- 
odromy data having two simple poles on the spectral plane. This monodromy data 
can be represented as a superposition of two one-pole terms and therefore, in ac- 
cordance with the previous section, it is naturally to expect that the corresponding 
solution will describe the superposition of fields of two Reissner - Nordstrom sources. 

Structure of the monodromy data 

Thus, we begin with the choice of the monodromy data functions in the form 

u{w) = — + — , v{w) = — + —. (20) 

w — hi w — h2 w — hi w — h2 

where Ui, U2, Vi, V2, hi and /12 are arbitrary complex constants. However, to obtain 
a static solution, we specify this data imposing the following constraints: 

ui = -Ml vi = -wT hi = /ii , , 

U2 = —U2 V2 = —V2 h2 = /l2 

Structure of the kernel and a new set of parameters 

In this case, the kernel W(r, C)/(C ~ ''') splits into the singular and regular parts: 



W(t,0 1 



c - T 



C-T {r-h)iC-h) 
whore the polynomials -P(C), R{t, C) defined as (everywhere below zq = 0): 

P(C) = hihl - 2hih2[{l - mi)/7.2 + (1 - tU2)hi]C + [{hi + h2f + 2/11/12 
- A{hiV2 + h2Vif - 2i{hfu2 + him) - Aihih2(ui + U2)]C 
+ [2/ii(-l + iui + 2m2 + AviV2 + 4w|) + 2/i2(-l + 2iui + iu2 

+ 4ViV2 + 4w2)](3 + [1 _ 2i{ui + U2) - i{vi + V2f] 

i?(T,C) =i?o(C)+i?i(C)r 

RoiO = hih2{hi + h2) 

- [hj + 3/ii/i2 + hj- ihi{hi + 2/12)^2 - ■i/i2(2/ii + /i2)wi] C 

+ [2(/ii -|- /i2)(l — iui — iu2 — 2W1W2) — ihiU2 — ih2Ui — Ahiv^ — Ah2vf\ Q"^ 

+ [-1 + 2i{ui + U2) + 4(ui + V2f] 

i?i(C) = —hih2 + {hi + /i2 - ihiU2 - ih2Ui)C + (-1 + iui + iu2)C 

(22) 

Assuming that the pole of the monodromy data is not located on the integration 
path L and therefore, hih2 0, we present P{w) in a factorized form 

P{w) = Po{w — 'Wi){w — wi){w — W2){w — W2) (23) 

The coefiicients of the polynomial P{w) always are real and therefore, its roots wi, 
wi and W2 , W2 are real or complex conjugated in pairs. We parameterize them as 

wi = zi + ai W2 = Z2 + (J2 hi = zi+mi 

, ,~ Z2 = zi+£ (24) 

Wi = Zi — (Ti W2 = Z2 — (72 h2 = Z2 + 1712 



February 2, 2008 3:44 



WSPC - Proceedings Trim Size: 9.75in x 6.5in main 



15 

where Zi and Z2 (and therefore, £) are real parameters while ai as well as (T2 can 
be real or pure imaginary. A comparison of coefficients of different powers of the 
spectral parameter in and in leads to a number of relations of the form 

2hih2{hi + h2~ ihiU2 - ih2Ui) PqIs: 

hl + hl+ Ahih2 - 2i{hi + h2){hiU2 + ft.2'"i) - 2i/?,i/i2(ui + U2) 

-A{hiV2+h2Vi)^ ^Pol2, (25) 

2{hi + h2) - 2i{hiU2 + h2Ui) - 2i{hi + /i2)(ui + "2) 
-8{hiV2 + h2Vi){vi + V2) = Poll, 

1 - 2i{ui + U2) - 4(wi + V2f = Po- 
where we have introduced the notations 

II = Wi + Wi + W2 + W2 

T2 = WiWi + W1W2 + W1W2 + W1W2 + W1W2 + ^2^2 
2^3 = Wi'W2W2 + W1W2W2 + W1W1W2 + WiWi'W2 
I4 = 'W1W1W2W2 

In terms of the parameters (|24p these functions possess the expressions: 

11 = 2{zi + Z2) la = 2^122(2:1 + Z2) - 2z2(tI - 2zict| 

12 ^ zl + 4ziZ2 + zl- cfI~ al X4 = {zl ~ af){z^ - (t|) 

The relations ((25|) allow to express the roots of P{w) (or the new parameters zi, 
Z2 and tJi, (72 in terms of ui, U2, wi, ^2 and hi, /12, however, it is more simple and 
convenient to use Zi, Z2 and ui, (T2 as new parameters and express Pq and some of 
the parameters ui, U2, wi, V2 and /ii, /i2 as functions of zi, Z2, <Ji, 02 and others. 
We present these expressions in the form 

_ 2ihihleie2PiP2 + 2ihmP^ + ihl{hi - h2){hlli - 2/1^X2 + 3/12X3 - 4J4) 

" 2(/li - /l2)3X4 ' 

2ihlh2eie2PiP2 - 2ihl^iPl + i/i^(/ii - h2){hlJi - 2/1^X2 + S/iiXg - 4X4) 



"2 



2(/li - /l2)3X 



4 



ih2eiPi ihie2P2 

^1 = -7:7; , , rr^ , V2 



2{hi~h2)^i 2{hi-h2)^i 

(26) 

where the parameters ei, 62 and Pi, P2 are defined as follows. Solving (|25|) leads to 
the expressions P(hi) and P(ft-2) which must be non-negative: 

P(/ii) EE (/ii - Wi){hi ~ Wi)(hi - W2){hi - W2) > . . 

P{h2) = {h2-Wi){h2-Wl){h2-W2){h2-W2)>Q ^ ' 

Assuming for the case of real roots of P{w) that they are numbered so that wi < 
wi < W2 < 1102, and that in any case, zi < Z2, we conclude from ([27]) that hi and 
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/i2 should be located outside the interval {wijWi) if wi (and therefore, wi) is real 
and outside the interval {w2,W2) if W2 (and therefore, •W2) is real. This means that 

{hi - wi){hi - wi) = ml - al > 0, {h2 - W2)(/i2 - W2) = ml - > 

and therefore, there exist the real parameters ei and 62 such that 

~2 ~2 2 ~2 ~2 2 /'oo\ 

and we allow to each of the parameters ei and 62 to be positive or negative in 
order to cover all opportunities to have in the expressions different signs in front of 
the square roots: ±^yP{hl) and ±^ P{h2)- The parameters Pi and P2 in (pS)) are 
defined as positive parameters which have the following explicit expressions: 



Pi = ^J{hl - W2){hi - W2) = yj{(.-mi )^ - cr. 



P2 = ^/{h2-wi){h2-wi) = y/{e + m2) 



2 _ ^2 

(29) 



The expressions (j26l) allow us to use instead uq, ui, vq, vi and hi, h2 the parameters 
zi, Z2, fhi, TO2, ei, 62 (with £ = Z2 — zi > 0) which can take arbitrary real values, 
provided Pi > and P| > 0, in view of (p9)) and (p8)) . Without loss of generality, 
we assume that if all roots are real, then cti > 0, (T2 > 0, and for the sources to be 
separated, we assume also that £ > ai + a2- 



Formal construction of the solution 

For solving of the integral equations ([8]) for the monodromy data ([20|l - ([2T|) . we 
substitute this data into the integral equation ([8]) and obtain the equations whose 
solution has the structure which is very similar to that for the one-pole case (fT4| : 

/<pf^'H\ / uw , . /Xa + Xiw + X2W^\ 

Vcp M (w)/ ^^""^ V Zo + Ziw + Z2W^ J 

where the coefficients Xk, Yk and Zk {k = 0, 1,2) are independent of the spectral 
parameter w and they are functions of coordinates and constant parameters. To 
find these coefficients explicitly, we have to substitute ([50)1 back into the integral 
equations, to calculate the corresponding integrals using the rules and solve the 
linear algebraic equations for Xk, Yk and Zk- Using this solution, we find explicitly 
(using again the rules (|16p ) all components ([6]) of the matrix R and calculate then 
pure algebraically the Ernst potentials and all metric and electromagnetic potential 
components in accordance with their general expressions ([5|). 



Weyl cylindrical and bipolar coordinates 

During the calculation of the contour integrals ([TS]) . the coordinate dependence of 
the solution arises from calculation of residues of the integrands. It is easy to see, 
that the residue at infinity gives rise to the terms which are polynomial functions of 
Weyl coordinates p and z, while the residues at finite poles — the roots of P(w) give 
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rise to the terms which are proportional to the values of the function X{^,r],w) at 
these roots: Ai = A(^, 77, wi), Ai = A(^, rj.wi), A2 = A(^, 77, W2) and A2 = A(^, 77, W2). 
To have deal with these functions, we introduce bipolar coordinates {xi,yi) and 
(2^2, J/2) related to the corresponding pairs of spherical- like coordinates (ri, 61) and 
{i'2,Q2) centered on the symmetry axis p — ed z — zi and z = Z2, so that 
xi — ri — mi, j/i — cos 01, X2 — r2 — rh2, 2/2 ~ cos6'2. The coordinates {xi, yi) and 
{x2,y2) are defined by the relations to Ai, Ai, A2 and A2 



^ xi - (Jiyi _ xi +aiyi J p = - crl^/l ~ 



Zl + (Tl Zl — CTl {z — zi+xiyi 

X2 - q-2j/2 r _ X2+ (T2y2 ( p ^ y^xf^~al^/i~^ 



A2 — ; A2 — 

Z2 + ^2 Z2 - 0-2 IZ = Z2+ X2y2 

while their relations to the Weyl coordinates given just above follow from the defi- 
nition of the function A(^, 77, w). 

Physical parameters of the solution 

The calculations described above lead to rather complicate form of the solution 
expressed in terms of functions of bipolar coordinates xi, X2, j/i, j/2, the formal 
mass and charge parameters ttt-i, 7712, ei, 62 and the parameters zi and Z2 which 
characterize the location of the sources on the symmetry axis. However, it is more 
convenient to describe the structure of this solution using certain combinations of 
these parameters, which, as it will be shown later, represent physical parameters - 
the masses of the sources tt^i, 7712, their charges ei, 62 calculated using the Gauss 
theorem and therefore, expressed in terms of the Komar-like integrals. These phys- 
ical parameters are determined by the expressions 

~ eie2(^ - mi -(- 7712) eiP2(^ - 777i + 7712) 

mi = 771-1 H -i^TT^ ^r-;^ ei — 



P1P2 + 6162 P1P2 + 6162 

eie2(£ ~ rhi + 7712) eiPi -I- e2P2 

7772 = 7712 p p , ~ ~ 62 = — ^ " 61 

^1^2 + 6162 i:-mi-|-m2 
where Pi, Pj were introduced in (j29p . The inverse relations are (7 is defined below): 

{£ + mi + m2)^ - Ko ^ [mi + m2f - P + Kq 

7771 — -, ^ m2 = ^ 

2{£ + mi+m2) 2{i + mi + m2) 

= - [(^ + 7772)2 -ml- Ko] + ei(ei + 62) 

^ = mi)2 - ml - Ko] + (62 + 27)(6i + 62) + ^T^^"!^^^" 

2 '■ ■' I + mi -t- m2 

Pi = i [(^ - -1)^ - + Ko] + (62 + 27)(61 + 62) - ^7;+";^^° 

2 I + mi + m2 

P| - ^ [(^ + ^2? - mf + Ko] + 61 (ei + 62) 
' {£ + m2)'^ - ml - Ko 



6162 = -PlP^ 



{£ + m2)'^ - ml + Ko , 
Kl = [{I + m2f - mlf + Aei{l + mi + 77t2)[62(^ - mi + m2) + 2eim2] 
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The parameters a\ and a-i also can be expressed in terms of mi, m-i, ei, €2 and I: 

i = ^r1 = "r1''?'^' where 7 ^ 7"' ' (31) 
crj = m2 — 62 m2 — 63 — 2e27, £ + mi + m2 

Using these parameters, we obtain more symmetric and simpler form of the solution. 
We introduce also, instead of ei and 62, two other parameters qi and q2 such that: 

m2qi-miq2 erf ^ mj + - qf, 

^= 'o ' 2 2,22 y^^' 

92=62+7, 1 0-3 ="^2+7 92- 

Finally, after some long calculations and consideration of various possible forms of 
our solution, its rather short form was found. In this form, all components of the 
solution are presented as functions of six parameters with only one constraint 



{mi, m2, gi, 92, £, 7} 



7^ = m2i7i - mi(j2 (33) 



which leaves only five parameters to be independent. 

5-parametric solution for interacting Reissner-Nordstrom sources 

For the monodromy data ((20|) . (|2T|) . the corresponding solution is static and its 
metric and electromagnetic potential in cylindrical Weyl coordinates take the forms 

At = Ap = = A,^ ^ 0, 

where iJ, / and $ are real functions of p and z. The calculations described in 
previous subsections lead to the following structure of the functions , $ and /: 

{V + gy ' ^ V + g' ^ {x\-oly\){xl-alyl) ^'"^ 

where 2?, t/, T are polynomial functions of bipolar coordinates with rather simple 
coefficients depending on the parameters of the solution: 



8\x\^ x\- a\y\ - a\y\ + 2(mim2 - 9192)2/12/2] 

g = mia;2 + m2a;i + 7(912/1 + 922/2) /^g-j 
+ 2(5 [mia;i + m2a;2 + 2/1(927 - "f^i^) + 2/2(917 + '^12^)] 

T = 91x2 + 922:1 + "fimiyi + m22/2) 

+ 2S[qiXi + q2X2 + 2/i('7i27 - 9i^) + 2/2("ii7 + 92^)] 

In psp , /o is an arbitrary constant which should be chosen so that / ^ 1 at spatial 
infinity and the parameter S in (j36p is determined by the expression: 

f - 1 , _ mim2 - 9192 , 

(1 + 25)2' ^2_^,2^^2^^2^^2 l-^U 
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In accordance with ([5^ and ([55]) . these equations give us an expHcit expression of 
the solution in terms of five free real parameters mi, TO2, ei, ei and i — zi — z\. 
We consider this solution as depending on six parameters mi, m2, gi, 92, ^ and 7 
restricted by the only one constraint (|33p and with the expressions ([52)1 for ct^, (t|: 

7^ = m2(7i - mi(72, (Ti = m^ + 7^ - , cr| = m2 + 7^ - g|. (38) 

This solution is asymptotically flat. As it will be explained in the next section, the 
parameters toi and TO2 characterize the individual masses of the sources and the 
charges of these sources are e\ =91+7 and 62 = 92 — 7 respectively, while the total 
mass m = mi + TO2 and the total charge e = ei + 62- 



Physical parameters of the sources 

For electrovacuum space-time which admits a time-like Killing vector field ^ = c?t we 
can write (following, for example,'^) the "dynamical" part of the Einstein-Maxwell 
equations in the Kinnersley-like self-dual form 

= 0, ^^r^ = (39) 
where the bivectors Ti^fc and Tn^ are defined as follows (7 = c = 1): 



J'ik — Fik + -£iklmF Fik — diAk — dkAi 



(40) 



and, due to ([39|) . these self-dual bivectors possess complex vector potentials, one of 
which determines the complex scalar function $ which enters the expressions (|40p : 



mk ^ d.Uk - ^kH^ T,k^^,'^k-^k<^^ $ = e*=$fc (41) 

The " dynamical" equations ([55)1 allow us to construct the " conserved" quantities 
- the additive integral values which characterize the sources and which can be 
calculated as the integrals over the spherical-like 2-surfaces surrounding different 
parts of the sources on the space-like hypersurfaces slicing the space-time region 
outside the sources. For stationary axisymmetric spacetime outside the field sources 
we consider the space-like hypersurfaces T,t : t = const where t is the Killing 
parameter. Contraction the equations ([55]) with the gradient dit, we obtain 

(3) (3) 

Vs{N,W^) = 0, Vs{N,r^) = (42) 

(3) , 

where 5 = 1,2,3; V is a covariant derivative with respect to a three-dimensional 
metric on the hypersurfaces St and Ni is a unit time-like future-directed nor- 
mal to these hypersurfaces. Integrating the above equations on Sf over a three- 
dimensional region between a sphere of a large radius Boo located in the asymptot- 
ically Minkowski region and by a closed surface Bs surrounding the sources, and 
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applying the Gauss theorem, we obtain 



(43) 

where Uk is a unit vector tangent to the hypersurface St and representing the 
outward normal to the corresponding boundary Bg or i?oo and d^a is the area 
element on these boundaries. At first, we calculate the integrals over B^o- 

For any stationary axisymmetric asymptotically flat electrovacuum solution of 
Einstein - Maxwell equations, the metric and electromagnetic potential of the form 
^ at spatial infinity admit the expansions (the NUT parameter, the total magnetic 
charge and the additive constants in At and are assumed to vanish): 

2m 1 , „ 2m „, 1 ^ 

= — + 0(^) f = '^ + — + 0{^) 

g,^ = ^sin2 + O(4) At = - + 0{-) 

g^^ = -r'^sm^e- 2mr sm^ e + 0{r°) A^ = 0(i) 

where p = r sinO, z = z^,-\-r cos 9 and the constants m, a and e mean the total mass, 
total angular momentum per unit mass and the total electric charge respectively. 
For the complex vector potentials Tii and $i introduced in ([4T|) for r ^ oo in the 
coordinates {t, p, z, if} we obtain the expansions 



2m 1 1 

Hi^i — — + 0{^), 0, 0, -2imcos6' + 0(-)} 

<i>,^{- + 0{\), 0, 0, iecos6' + 0(-)} 
The components of the vector N'^ in these coordinates take the form 



(44) 



^' = {i + 7 + o(^), 0, 0, ^ + o(-l)} 

For a sphere of a large radius, the spatial unit normal vector in the leading order is 
n^di ~ d/dr and in the limit r ^ oo for the integrals (|43p over Boo we obtain 

j {HikN'n'')d^a = -Snm, J {T^kN'n^)d^a = ^-kc 

Boo Boo 

where m and e are the total mass and charge of the field configuration. This allows 
us to express (in accordance with ((43|) ) the total mass and charge of a system of 
sources in terms of the integrals over the surface Bs surrounding the sources: 

m = j {n^kN'n^)d^a, e = ^ j {J',kN'n^)d^a (45) 



Bs 
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+ + 

Now we transform the integrands in (|45)l using the self-duality of Ti and T: 



B^, B, 



e'l^^Nkui = k'p - yr 



(46) 



Here we introduced two spatial unit vectors fc* and which are orthogonal to each 
other and tangent to 2-surface Bg. These vectors determine the orientation of the 
2-surface Bg so that the spatial basis {n', y , l'^} is positively oriented. Choosing 
the 2-surface Bg to be axially symmetric, we specify the choice of the vectors fc' 
and P in the tangent space of Bg so that P would be the rotational Killing vector 
d/d^ and k^ is tangent to a curve C on Bg with ip = const: 



dp dz 



1 



{0, 0, 0, 1} 



where the parameter £ is the length on the curve C and the direction on C is chosen 

so that it goes from some point on the positive part of the z-axis to some point 

on its negative part (provided the point p = z = is located somewhere inside 

Bg). The element of the area on Bg is y/^g^d£dip. After the integration over ip 

we observe that the integrands of the remaining contour integrals over C are the 

-I- + 

differentials of the (^-components of the potentials of self-dual bivectors Ti. and JF: 



(47) 



where (. . .)+ means the value of a potential at the beginning of C, i.e. on the positive 
part of the z-axis, and (. . .)_ ~ its value on the negative part of this axis. 

It is important to note that Ti.^ and are constant along the regular parts 
of the symmetry axis outside the sources and therefore, the integrals (|47l) do not 
change if we deform the surface of integration Bg surrounding the sources in the 
domain of regularity of the solution. 

Now we consider the case, in which the source of the field consists of two parts 
(black holes, naked singularities or extended bodies) separated by some segment 
of the symmetry axis. In this case, we can deform the surface Bg surrounding the 
sources into two spheres each surrounding one of these parts of the sources and 
a thin tube surrounding the intermediate part of the axis between the sources. If 
7i(^ and are constant on the intermediate part of the axis, the integral over 
the thin tube vanish and we obtain that the integrals for the total mass and the 
total charge are expressed as the sums of the integrals of the same type over the 
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spheres surrounding each of the sources. In particular, for the sources consisting of 
two parts, as it is in our solution for two Reissner - Nordstrom sources, we obtain 

TO — TOi + TO2, e — ei + 62 

where mi, m2 and ei, 62 can be interpreted respectively as the masses (energies) 
and charges of the corresponding parts of the source. To conclude this section, we 
mention that rather tedious calculations show that the parameters toi, TO2 and ei, 
62 in our electrostatic solution presented in the previous section and describing the 
field of two Reissner - Nordstrom sources, coincide with the values of the integrals 
([T7|) calculated for each part of the sources and therefore, these parameters can be 
interpreted as the additive masses and charges of these interacting sources. 

Equilibrium of two Reissner - Nordstrom sources 

The 5-parametric solution ([34| - ([38|) presented above is asymptotically flat and 
the space-time geometry is regular far enough from the sources. However, for an 
arbitrary choice of parameters toi, TO2, qi, 92 and £, this solution may have some 
physically not reasonable singularities (" struts" ) on the part of the axis of symmetry 
between the sources. Like at the vertex of the cone, at these points, the local Eu- 
cledean properties of space-time may be violated so that on the space-time sections 
t — const, z — const the ratio of the length £ of a small circle, surrounding the 
axis of symmetry p — and contracting to the point z of this axis, to its radius TZ 
multiplied by 27r is not equal to a unit. The constraint imposed on the parameters 
which provide the absence of such conical singularities plays the role of equilibrium 
condition because it allows to select a physically acceptable solution in which these 
singularities are absent and the sources are in the equilibrium because of the balance 
between their gravitational and electromagnetic interactions. 

As it follows from elementary considerations, the limiting value of the square of 
the ratio {27tTZ) /£ for 7?. — > is equal to the value of the product fH at the point 
z of the axis. However, the field equations imply that on the regular parts of the 
axis of symmetry the product f{p = 0, z)H{p = 0, z) does not depend on z and 
therefore, it is a constant. However, the values of this constant can be different on 
different parts of the axis of symmetry separated by the sources, and the condition 
fH = 1 of the absence of conical singularities should be satisfied at each of these 
parts of the axis of symmetry. 

In our solution with two separated sources, there exists three parts of the axis of 
symmetry where the condition of the absence of conical points should be considered. 
There are two semi-infinite parts of the axis outside the sources. On the negative part 
L_, we have —00 < 2 < zi — cti (if the first source is a black hole) or —00 < z < zi 
(if the first source is a naked singularity) and on this part xi ^ zi — z, X2 = Z2 — z, 
2/1 = —1, y2 — —1. On the positive part L+ of the axis, Z2 + 0^2 < z < 00 {if 
the second source is a black hole) or Z2 < z < 00 (if the second source is a naked 
singularity) and we have there xi = z — zi, a;2 = z — 22, J/i = 1, J/2 = 1- It is easy to 
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see that on these parts of the axis the condition fH = 1 is satisfied for any choice 
of free parameters of our solution. However, it is not the case on the part Lq of 
the axis between the sources. On Lq we have zi + Re{ai) < z < 22 ~ Re{a2) and 
xi = z — zi, X2 = zi + i — z, Hi = 1, ?/2 = —1 and the product fH also takes there 
some constant value whose equality to a unit give us the equilibrium condition: 

TOir7i2 — qiq2 = (48) 

It is interesting to note that this equilibrium condition looks just like the Newtonian 
condition of equilibrium of two charged point-like masses, but in the case of General 
Relativity this condition relates the masses of the Reissner-Nordstrom sources not 
with their charges ei and 62, but with the parameters qi and q2 whose expressions 
in terms of masses mi, m2 and charges ei, 62 depend also on the z-distance £ 
separating the sources. 

The equilibrium condition (|48| allows also to simplify the solution p4|) - (|38p 
and leads to the 4-parametric family which describes the superposed field of two 
Reissner - Nordstrom sources in equilibrium. This solution has been presented in 
our previous short paper. ^ Here we persent it in a bit different form using mi, TO2, 
9i, 927 7 as the basic set of parameters, such that the first four of them should 
satisfy the equilibrium condition (|48p . For this solution the metric functions have 
the same expressions (|35|) where the polynomials 2?, Q and T have the expressions 



mim2 = (7192 



V = X1X2 - 7^2/12/2 
g = mia;2 + m2Xi + 7(911/1 + 922/2) 
T 91X2 + 922^1 + j{miyi + 17121/2) 

and /o = 1. The other physical parameters - the z-distance i separating the sources 
and their charges are determined by the expressions 

£ = Z2 - 21 = (m29i - mi(j2)/7, 61=91+7, 62=92^7 

The independent parameters mi, m2, 91, 92, 7 of this solution should be chosen 
so that the sources would be separated actually by some positive distance, i.e. 
i > Re{(Ti)+Re{(J2). 

This concludes our description in this paper of some fragments of the monodromy 
transform approach and of the procedure for solution of the corresponding singular 
integral equations for stationary axisymmetric electrovacuum fields with simple ra- 
tional monodromy data, which leads to the construction of the 5-parametric family 
of solutions for the field of two interacting Reissner-Nordstrom sources, and of the 
derivation from this solution of the 4-parametric family of solutions for the fields 
of equilibrium configurations of these sources. Some properties of these equilibrium 
field configurations have been discussed in,^ however, a more detail analysis of phys- 
ical and geometrical properties of these configurations, such as the structure of the 
superposed fields, influence on the geometry of horizons and on the space-time ge- 
ometry inside the horizon of the external gravitational and electromagnetic fields 
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created by another source, tidal influence of these fields on the structure of naked 
singularities, stability of equilibrium and probably, some others, are expected to be 
the subject of our next publications. 
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